Astronomy & Astrophysics manuscript no. GonzalezPerez' Metcalfe 


© ESQ 2008 


Febraary 15, 2008 





Testing the evolution of the DB white dwarf GD 358: First results 
of a new approach using asteroseismology 

Jose Miguel Gonzalez Perez ^'^ and Travis S. Metcalfe^ 



Instituto de Astrofisica de Canarias, 38200 La Laguna, Tenerife, Spain 
^ 2 GRANTECAN S.A. (CALP), 38712 Brefia Baja, La Palma, Spain 

^ High Altitude Observatory, National Center for Atmospheric Research, P.O. Box 3000, Boulder, Colorado, USA 

o 

Received (date) / Accepted (date) 
^ ABSTRACT 

Aims. We present a new method that investigates the evolutionary history of the pulsating DB white dwarf GD 358 using asteroseismology. 
^ This is done considering the internal C/0 profile, which describes the relative abundances of carbon and oxygen from the core of the star to its 
surface. Different evolutionary channels lead to the generation of different C/O profiles, and these affect the pulsation periods. 
1—1 Methods. We used C/O profiles associated with white dwarfs that evolved through binary evolution channels where the progenitor experienced 
one or two episodes of mass loss during one or two common envelope (CE) phases, and two profiles from single star evolution. We computed 
CIh models using these different profiles and used a genetic algorithm (GA) to optimize the search in the parameter space in order to find the 
Q best-fit to the observed pulsation periods. We used three-parameter models, adjusting the stellar mass (M*), the effective temperature (Teff), 
5_j and the helium mass of the external layer (Mh,,). 

Results. Our results suggest that binary evolution profiles may provide a better match to the pulsation periods of GD 358. The best-fit to the 
Ci observations is obtained using a profile related to an evolutionary history where two episodes of mass loss happen during two CE phases, 
'— ' the first during the RGB (Red Giant Branch) stage. The values obtained are r^//=24 300 K , M,=0.585 M© , and log(MHjM*)=-5.66. The 
_ ^ best-fit model has a mass close to the mean mass for DB white dwarfs found in various works, and a temperature consistent with UV spectra 
^ obtained with the lUE satellite. 
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culated accurately by evolving simpler, less detailed models. 
The DAV stars are generally modeled with a core of carbon 



potential amount of information encoded in their pulsation fre 
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lium layer surrounding the carbon/oxygen core. Robinson et al. 
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Whole Earth Telescope (Nather et al. 1990). 

There are presently three known classes of pulsating white 
dwarfs. The hottest class is the planetary nebula nucleus vari- 
ables, which have atmospheres of ionized helium and are also 
called DOVs. These objects show complicated and variable 
temporal spectra (Gonzalez Perez et al. 2006), and require de- 
tailed calculations that evolve a main sequence stellar model 



The chemical evolution of the DB white dwarfs atmo- 
spheres cannot be satisfactorily explained to date (see Shipman 
1997 for a review). In particular, the presence of the DB gap, 
which is the absence of DB white dwarfs between tempera- 
tures of -30 000 - 45 000 K (e.g. Liebert 1986), is pooriy un- 
derstood. This DB gap persists in the latest results from the 
Sloan Digital Sky Survey (Kleinman et al. 2004). Eisenstein et 
Send offprint requests to: J. M. Gonzalez Perez, e-mail: al. (2006) found several DBs in the gap, but still a decrease in 
j gperez®iac . es the relative number. One possibility is that DBs might turn into 
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DAs (white dwarfs with pure H atmospheres) through diffusion 
of small traces of H to the surface of the star, concealing the He 
atmosphere. However, Castanheira et al. (2006) showed that at- 
mospheric contamination with H is not directly proportional to 
Teff for DB stars, in conflict with what would be expected un- 
der the diffusion scenario. 

It is also not clear whether DBs are mostly produced by 
single-star evolution or whether a significant fraction of them 
originate from binary progenitors. Nather et al. (1981) pointed 
out that interacting binary white dwarfs (IBWDs) should pro- 
duce DBs at the end of their evolution. AM CVn, one of the 
best studied IBWDs, has an T^ff ~ 25000 K. This tempera- 
ture is inside the interval where DBVs are found. If AM CVn 
ends mass transfer phase at this temperature, there will be a 
single DB left, quite possibly a DB V. The examination of these 
problems can be aided by asteroseismology. Compared with 
other classes of pulsating stars, the DBVs are promising can- 
didates for asteroseismology: their mode spectra are relatively 
rich, and pulsation theory for these stars is quite advanced and 
well tested. 

GD 358, also called V777 Hercuhs, is the prototype of the 
DBVs and is one of the classical examples of the appUcation 
of asteroseismological methods (e.g. Winget et al. 1994; Vuille 
et al. 2000; Kepler et al. 2003). GD 358 was the first pulsat- 
ing star detected based on a theoretical prediction (Winget et 
al. 1985). When the Whole Earth Telescope observed GD 358 
in 1990, 181 periodicities were detected, but only modes from 
radial order A: = 8 to 18 were identified, most of them showing 
triplets, consistent with the degree I - 1 identification (Winget 
et al. 1994). Detecting as many modes as possible is important, 
as each independent mode detected yields an independent con- 
straint on the star's structure. The authors derived the total stel- 
lar mass as 0.61 + 0.03 M©, the mass of the outer hehum enve- 
lope as (2.0± 1.0)x 10-^ M*, the luminosity as 0.050±0.012 Lq 
and, deriving a temperature and bolometric correction, the dis- 
tance as 42 ± 3 pc. 

Asteroseismological analysis of GD 358 has been a suc- 
cessful tool to test pulsation theory. For instance, the cen- 
tral oxygen abundance of the models has been shown to have 
a measurable efl'ect on the pulsation frequencies (Metcalfe 
et al. 2001 -hereafter Ml-), and also the possible presence 
of a ^He/'*He transition zone caused by chemical diffusion 
(Montgomery et al. 2001), although more recent observations 
have ruled out this possibility (Wolff et al. 2002). 

2. Pulsations and C/0 profiles. 

The interior structure of the white dwarfs contains a record of 
the physical processes that operate during the later stages in the 
lives of most stars. The distribution of the C/O ratio depends 
on the evolution of the star. Stellar evolution theory predicts 
that a large number of white dwarfs are post Asymptotic Giant 
Brach (AGB) stars. Therefore, they should consist of the pri- 
mary ashes of He-burning, i.e. carbon and oxygen. The central 
ratio C/O depends on the interplay between the triple-a and the 
'^C(q', y)'*(9 reactions during the core helium burning phase. 
In the first part of the burning the most efficient process is the 
triple-or, which produces ^^C. When the central abundance of 



heUum decreases below about 0.15 Mi,, the production of °0 
via the ^^C(a, 7)'^0 becomes the dominant process. 

While they are still embedded in the cores of red giant mod- 
els, the internal chemical profiles of white dwarfs models show 
a relatively constant C/O ratio near the center, with a size deter- 
mined by the extent of the helium-burning convective region. 
Further out the oxygen mass fraction decreases as a result of 
the helium-burning shell moving towards the surface of the red 
giant model while gravitational contraction causes the temper- 
ature and density to rise. This increases the efficiency of the 
triple-a reaction, producing more carbon relative to oxygen. 
The central oxygen mass fraction is lower in higher mass white 
dwarfs models. The rate of the triple-a reaction increases faster 
at higher densities than does the '^C(a, y)'^^ reaction. As a 
consequence, more helium is used up in the production of car- 
bon, and relatively less is available to produce oxygen in higher 
mass models. 

Metcalfe et al. (2000 -hereafter M0-) used the observed pe- 
riods of GD 358 from Winget et al. (1994) to search for the op- 
timum theoretical model with static diffusion envelopes. This 
search was based in the application of a genetic-algorithm- 
based optimization (GA) to white dwarf pulsation models, 
which performs a global search to provide objective global 
best-fit models for the observed pulsation frequencies. They 
used a three-parameter model including Mi,, Tgff, and Mhc, 
and six different combinations of core composition and C/O 
transition profiles to search for the best-fit model. Their results 
showed that both the central oxygen abundance and the shape 
of the C/O profile affect the pulsation pattern. 

Ml modified the code to include as free parameters any 
central oxygen mass fraction (Xq) between and 1 with res- 
olution 1% and a fractional mass parameter {q). The GA fit- 
ting process explored different chemical profiles built in the 
following way: the value Xq was fixed to its central value out 
to q that varied between 0.1 and 0.85 with resolution 0.75%. 
From this point Xq decreased linearly in mass to zero oxy- 
gen at the 95% mass point. The results showed that signifi- 
cant improvement to the fit was possible by including Xq and q 
as free parameters, confirming that the observed pulsations re- 
ally contain information about the hidden interiors of the stars. 
The best fit obtained was: T^fj = 226QQK; Mi, = O.65OM0; 
\og{MHelMi,) = -2.74; Xq = 0.84; q = 0.49. 

Therefore, the C/O profiles may give us an opportunity to 
explore the previous evolution of the DBV white dwarfs. C/O 
profiles obtained from single evolutionary models are differ- 
ent than those obtained from evolutionary models of close bi- 
nary systems (e.g. Iben & Tutukov 1985 -hereafter 1T85-; Iben 
1986). The shapes of the C/O profiles may be affected due to 
interaction with close companions, i.e., if CE phases happen 
during the RGB and/or AGB stages. 

3. Close binary evolution. 

Peculiar white dwarf stars may be created through binary 
evolutionary channels. We have concentrated on evolution- 
ary channels that produce white dwarfs without hydrogen. 
Hydrogen may be removed of the system due to previous CE 
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stages. The resulting white dwarfs may have a surface rich in 
helium and may contribute to the DB population. 

IT85 analyzed the possibility to obtain such white dwarfs 
by evolving components of close binary systems where the 
component masses were in the range 2 to 12 Mq. They as- 
sumed models where the more massive component filled its 
Roche lobe during the RGB, after the central hydrogen exhaus- 
tion, but before or near the onset of core helium burning. The 
binary components experience one or two mass loss processes 
due to CE stages, depending upon the initial main-sequence 
mass. 

Components with masses in the range 2.3 to 4.8 Mq expe- 
rience one episode of mass loss and become degenerate dwarfs 
of mass in the range 0.32 to about 0.7 Mq. The remnants ex- 
perience several small hydrogen shell flashes and one large fi- 
nal flash which, for the more massive degenerates dwarfs, re- 
move almost all of the hydrogen-rich layer. In these cases the 
total mass of the hydrogen layer is reduced to about IQ~^ Mq. 
This evolutionary channel may produce remnants with masses 
as small as 0.4 Mq and can convert 70% of the initial helium 
into carbon/oxygen. This represent a substantial difference 
from evolutionary channels of single low-mass stars, which 
do not ignite and burn helium until the mass of the electron- 
degenerate cores reaches ~0.5 Mq. Han et al.(2000) have per- 
formed close binary evolutionary calculations and found that a 
carbon/oxygen white dwarf with mass as low as 0.33 Mq may 
be formed from stable Roche overflow if the initial mass of the 
primary is close to 2.5 Mq. 

IT85 found that binary components of initial mass in the 
range 4.8 to 10.3 Mq lose the hydrogen-rich matter in two 
episodes of Roche-lobe overflow and become C/O degenerate 
dwarfs with masses in the range 0.7 to 1.08 Mq. This chan- 
nel leads to the production of DB white dwarfs. For the more 
massive remnants (over 0.8 Mq), a substantial part of the he- 
lium layer is also removed from the system during the second 
Roche-lobe overflow, so the degenerate dwarf would have a 
thin helium outer layer. This layer is expected to be thicker for 
remnants below 0.8 Mq . IT85 estimated a formation rate of 
~ lO^^yr ' of non-DA degenerate dwarfs created for these bi- 
nary scenarios and concluded that they may produce 15%-20% 
of all hot non-DA white dwarfs. 

4. Computation of models. 

4.1. DB white dwarfs models. 

To find the theoretical model of a white dwarf, we start with a 
static model of a pre-white dwarf and allow it to evolve quasi- 
statically until it reaches the desired temperature. We then cal- 
culate the adiabatic non-radial oscillation frequencies for the 
output model. It is not necessary to go through detailed calcu- 
lations that evolve a main sequence star all the way to its pre- 
white dwarf phase. For the DB and DA white dwarfs, it is suf- 
ficient to start with a hot polytrope of order 2/3 (i.e. P oc p^^^). 
The cooUng tracks of these polytropes converge with those of 
the pre-white dwarfs models above the temperatures at which 
DB and DA white dwarfs are observed to be pulsationally un- 
stable (Wood 1990). 



We have used the White Dwarf Evolution Code (WDEC) 
to evolve a start model to a specific temperature. This code is 
described in detail by Lamb & Van Hom (1975) and by Wood 
(1990). It was originally written by Martin Schwarzschild, and 
has subsequently been updated and modified by many oth- 
ers including: Kutter & Savedofl" (1969), Lamb & Van Horn 
(1975), Winget (1981), Kawaler (1986), Wood (1990), Bradley 
(1993), and Montgomery (1998). The pulsation frequencies 
of the output models are determined using the adiabatic non- 
radial oscillation (ANRO) code described by Kawaler (1986). 
Practical modifications to these programs were done by MO, 
primarily to allow models to be calculated without any inter- 
vention by the user. 

The equation of state for H/He comes from Saumon et 
al. (1995). We use the OPAL opacity tables from Iglesias & 
Rogers (1996), neutrino rates from Itoh et al. (1996), and the 
ML2 mixing-length prescription of Bohm & CassinelU (1971). 
We have fixed the ratio of mixing-length/pressure scale height 
to 1.25, as recommended by Beauchamp et al. (1995). The 
evolution calculations for the core are fully self-consistent, 
but the envelope is treated separately. The core and envelope 
are stitched together and the envelope is adjusted to match 
the boundary conditions at the interface. Adjusting the he- 
hum layer mass involves stitching an envelope with the desired 
thickness onto the core before starting the evolution. Because 
this is done while the model is still very hot, there is plenty of 
time to reach equilibrium before the model approaches the final 
temperature. 

4.2. The genetic algorithm and Darwin. 

The global search for the optimal model parameters to fit the in- 
dependent pulsations periods of GD 358 is done using the GA 
optimization method developed by MO and Ml. This method 
improves the objective search in the parameter space over tra- 
ditional procedures like iterative methods starting from a first 
guess. This standard approach has a potential problem: the ini- 
tial set of parameters is usually determined subjectively. This, 
combined with a local search in the parameter space, may pro- 
duce a local good fit that is not the global best fit. Restrictions 
by the GA method regarding the range of the parameters space 
are imposed only by observational constraints and the physics 
of the model. Therefore, GA provides a relatively efficient way 
of searching globally for the best-fit model. 

The first step of the GA is to fill up the parameter space with 
trials consisting of randomly chosen values for each parameter, 
within a range based on the physics that the parameters are sup- 
posed to describe. The model is evaluated for each trial, and the 
result is compared with the observational data. This compari- 
son assigns a fitness to the trial inversely proportional to the 
variance. The new generation of trials is then created taking 
into account the previous fitness. The fittest models survive to 
the next generation and the parameter space is then sampled 
more around these values. Operators emulating reproduction 
and mutation are applied to produce the new generation of tri- 
als. The evolution continues for a specified number of gener- 
ations, chosen to maximize the efficiency of the method. We 
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used a population size of 128 trials and allowed the GA to run 
for 200 generations. We have performed a total of 10 GA runs 
for each core composition to reduce the chances of not finding 
the best answer to less than about 3 in 10 000. 

The GA code has been implemented in a network of 64 PCs 
running Linux (Metcalfe & Nather 1999 and 2000). The white 
dwarf code runs on this metacomputer by using the Parallel 
Virtual Machine (PVM) software (Geist et al. 1994), which al- 
lows a collection of networked computers to cooperate on a 
problem as if they were a single multiprocessor parallel ma- 
chine. The central computer in the network is called Darwin 
and runs the parallel version of the genetic algorithm (PIKAIA) 
which is the master program. After creating a new generation, 
Darwin distributes to the slave computers an array with the 
models to check. The slave computers perform the following 
actions: evolve a white dwarf model to the specified tempera- 
ture, determine the pulsation periods of the model, and com- 
pare observations with the calculated pulsation periods. Then 
it sends the result back to Darwin which will include it in the 
GA code. The slave computer is then ready for a new trial so 
Darwin can send a new trial from the generation. 

4.3. Parameter space. 

We use a three-parameter model including M*, T^j j, and Mh^- 
The GA search space is defined as follows: the masses are 
confined between 0.45 Mq and 0.95 Mq. Although Kepler et 
al. (2007) reported a substantial number of white dwarfs over 
this upper limit, Beauchamp et al. (1999) and Castanheira et 
al. (2006) found that all known DBVs appear to have masses 
within this range. The temperature search includes values be- 
tween 20 000 K and 30 000 K and is based in the temperature 
determination for the known DB instability strip, that is, the 
strip in the T^jj versus log(g) diagram where the DB white 
dwarfs are pulsationally unstable (Beauchamp et al. 1999). 
These authors, depending on various assumptions, place the 
red edge as low as 21 800 K, and the blue edge as high as 
27 800 K. Finally, the search interval for the mass of the atmo- 
spheric helium layer has a lower limit 10"*^ M^, and an upper 
limit 10"^ Mi,. Masses must not be greater than this because 
the pressure on the overlying material will then theoretically 
initiate helium burning at the base of the envelope. At the other 
extreme, none of the models pulsate for helium masses less 
than 10"** over the entire temperature range we consider 
(Bradley & Winget 1994). 

A two-digit decimal is used for encoding the different pa- 
rameters. This results in a temperature resolution of 100 K, a 
mass resolution of 0.005 Mq, and a resolution for the helium 
layer thickness of 0.05 dex. 

4.4. Computing profiles. 

We have used the C/O profiles presented in IT85 and Iben 
(1986). We had to modify the original profiles to be able to in- 
corporate them into our code. The axes have been transformed 
to Xq (mass fraction of oxygen) versus Mr/M-t. (enclosed mass 
fraction at a given radius). A C/O profile is computed using 
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Fig. 1. Profiles of evolutionary channels affected by 2 epochs 
of mass loss during 2 CE phases, the first one during the 
RGB stage. The upper panel shows the profile obtained for a 
0.752 Mq remnant and the lower for a 0.894 Mq remnant. 
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Fig. 2. Profiles of an evolutionary channel affected by 1 CE, 
and 2 evolutionary channels of single stars. The upper panel 
shows the profile obtained with a evolutionary history affected 
by 1 process of mass loss due to 1 CE during the RGB stage 
for a 0.523 Mq remnant. The middle panel shows the pro- 
file obtained, using a single evolutionary model, for a white 
dwarf model of mass 0.599 Mq that evolved from a horizontal- 
branch star of initial mass 7 Mq. The lower panel shows an- 
other profile related to a single evolutionary scenario: a white 
dwarf model of mass 0.965 Mq that evolved from a single AGB 
model star of mass 7 Mq 
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a series of pairs of points (Xc, MrlMi,), where Xc - I - Xq. 
These pairs define the profile shape in the model. We have used 
between 6 and 8 points. This procedure simplifies the original 
shape of the profile but reproduces it well. One limitation for 
computing the profiles is that Xq must be set to zero outside 
the mass fraction (MrlM)-Q.95. However, some of the original 
profiles present some amount of oxygen outside this fractional 
mass. Another simplification is that our models consider car- 
bon and oxygen as the only nuclear species in the interiors of 
the white dwarfs. IT85 evolutionary models present traces of 
other elements like '^^Ne and ^^Mg. Figures [T] and |2] show the 
profiles we tested. 

Our study is mainly based on IT85 because, until now, it 
is the only work presenting several C/O profiles for different 
evolutionary channels of close binary systems. In addition, the 
profiles are obtained using the same physics, which improve 
the possibilities for a quantitative comparison. Other profiles 
for single evolutionary scenarios are available (e.g. Salaris et 
al. 1997). However, the physics included in these evolutionary 
models are slightly different, e.g. how overshooting is treated, 
nuclear reactions rates, etc. Models with different physics may 
produce different C/O profile shapes. 

IT85 evolutionary models start with an initial pair of stars 
with certain masses and separation and are evolved to produce 
a white dwarf with a certain C/O profile and mass, mainly with- 
out hydrogen on the surface. Therefore, each C/O profile tested 
is already related to a certain mass. In our code, M* is a free 
parameter. A self-consistent result for a given profile would re- 
quire a best-fit model with a value for M^, not too far from 
the final masses in IT85, though profiles for different masses 
should be qualitatively similar. 

Table [T]presents the properties of the different C/O profiles 
we tested. The second column indicates how many episodes of 
mass loss have happened due to CE phases in the previous evo- 
lution. RGB indicates that the first CE happened during the Red 
Giant Brach stage. SE indicates single evolution. HB indicates 
a white dwarf remnant coming from a horizontal branch model, 
and AGB from an Asymptotic Giant Branch model. We have 
included two profiles associated with a degenerate white dwarf 
that has experienced 2 epochs of mass loss in CE episodes, the 
first one during the RGB phase: profile 1 is from the 0.752 Mq 
remnant of a model of initial mass 5 Mq (see upper panel of 
Figure [TJ, while profile 2 comes from a 0.894 Mq remnant of 
a model with initial mass 6.95 Mq (see lower panel of Fig. [TJ. 
Profile 3 corresponds to the 0.523 Mq remnant of a model of 
initial mass 4 Mq that has experienced only one episode of 
mass loss due to a CE during the RGB phase. This model re- 
tains some hydrogen in the surface layer, ~2.6 x 10""* (see 
profile in the upper panel of Fig.|2]). 

We have completed our comparative analysis by checking 
the best-fit models for 2 profiles obtained from single evolu- 
tionary scenarios. Profile 4 is related to a white dwarf model of 
mass 0.599 Mq that comes from the evolution of a horizontal- 
branch model of initial mass 7 Mq (Iben 1982). This is pre- 
sented in the middle panel of Fig. |2]i. Profile 5 is the remnant 
0.965 Mq core of a single AGB model star of mass 7 Mq (Iben 
1976) (see lower panel of Fig.|2]i. This last model presents some 



inconsistency since our parameter space for the mass is limited 
to masses below 0.95 Mq. 



5. Results. 

Table |2]presents the best fits obtained for each profile. The best 
fits are obtained for the profiles 1 and 2, which involve evo- 
lutionary channels where the star experiences two CE phases. 
Profile 1 is associated in the IT85 evolutionary model to a white 
dwarf with mass 0.752 Mq, and profile 2 to a more massive 
white dwarf (0.894 Mq). The asteroseismological masses ob- 
tained are in both cases smaller: 0.555 and 0.585 Mq. Both fits 
have almost the same residuals, 2.13 and 2.12 s, but the values 
obtained for the 3 parameters, T^//, M*, and log(M//(./M*), are 
slightly different. This confirms the dependence of the models 
on the intrinsic shape of the profile. 

The best fit is obtained using profile 2. The parame- 
ters obtained are T^ff^lA 300 K , M*=0.585 Mq , and 
log(MH(./M*)=-5.66. The best-fit model has a mass close to the 
mean mass for DB white dwarfs, and a temperature consistent 
with UV spectra obtained with the lUE satellite (Castanheira et 
al. 2005). 

Figure |3] shows the general characteristics of the three- 
dimensional parameter space: M*, T^jj, and log(M//(./M*). All 
the combinations of parameters found by the GA algorithm for 
our best fit model with residuals smaller than 3.3 s are pre- 
sented. The two panels are orthogonal projections of the search 
space, so each point in the left panel corresponds to a point in 
the right panel. The fit is displayed using the following color 
scale: a red point if the residuals (r.m.s.) are smaller than 2.2 
s; yellow: 2.2 < r.m.s. < 2.4; green: 2.4 < r.m.s. < 2.7; cyan: 
2.7 < r.m.s. < 3.0; blue: 3.0 < r.m.s. < 3.3. The figure shows 
the presence of more than one region that yields a good match 
to the observations. Two families of good fits, for thick and thin 
helium layers, are displayed, but the amount and quality of the 
fits found for thin layers are larger. The distribution of the best 
models form stripes in the search spaces panels. This is due to 
the parameter correlations in both projections. The correlation 
between M* and log(M//(,/M*) is described by Brassard et al. 
(1992), who showed that a thinner helium layer can compen- 
sate for an overestimate of the mass. 

The correlation between and Teff is related to the 
Brunt-Vaisala frequency, which reflects the difference between 
the actual an the adiabatic density gradients. The pulsation 
periods of a white dwarf model reflect the average of the 
Brunt-Vaisala frequency through the star. If the temperature 
decreases, the matter becomes more degenerate, so the Brunt- 
Vaisala frequency tends to be zero in much of the star, leading 
to lower pulsation frequencies. So an overestimate of the mass 
can compensate for an underestimate of the temperature. 

Figure |4] presents the observed and calculated periods for 
the best model plotted against the forward period spacing {6P = 
Pk+\ - Pk)- Note that the GA only fits the periods of the pul- 
sations modes, and the agreement between the deviations from 
the mean period spacing is a reflection of the overall quality of 
the match. 
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Table 1. Characteristics of the evolutionary channels of the profiles we tested. The second column indicates how many episodes 
of mass loss have happened due to CE phases in the previous evolution. RGB indicates that the first happened during the RGB 
stage. SE indicates single evolution. HB indicates that the white dwarf evolved from a Horizontal Branch model, and AGB from 
an Asymptotic Giant Branch model. 



Profile id 


Evolution 


Mass WD (Mq) 


Hydrogen 


1 


2 CE, RGB 


0.752 


no 


2 


2 CE, RGB 


0.894 


no 


3 


1 CE, RGB 


0.523 


~ 2.6 X 10-" 


4 


SE, HB 


0.599 


yes 


5 


SE, AGB 


0.965 


yes 



Table 2. Best models obtained for the different evolutionary channels we tested. The fits labeled with 'a' are not compatible with 
the observed parallax of GD 358, which constrains the luminosity. 



Profile id 


T.ff (K) 


M/Mq 


log(M„,/MJ 


r.m.s. (s) 


1 


25 000 


0.555 


-5.60 


2.13 


2 


24 300 


0.585 


-5.66 


2.12 


3 


23 700 


0.605 


-6.03 


2.19 


4 


25 100 


0.825 


-3.96 


2.23° 


5 


26 000 


0.560 


-5.39 


2.21 



0.9 



O.B 



\ 0.7 h 



0.6 - 



0.5 
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Temperature (K) 
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Fig. 3. Front and side views of the search space for the best model (profile 2). Square points mark the location of models that 
yields a reasonable match to the periods observed for GD 358. The color of a point indicates the relative quality of the match (see 
text for details). 



6. Discussion. 

We have presented a new method to investigate the evolution- 
ary history of GD 358, which uses the dependence of the pul- 
sation frequencies on the internal C/O profile. Different evo- 
lutionary channels lead to the generation of different profiles 
describing how the C/O ratio changes from the core of the star 
to its surface. We have tested and compared profiles associated 



to binary evolution where the object experiences one or two 
episodes of mass loss during one or two CE phases. In all the 
cases the first CE happens during the RGB phase. We have also 
analyzed two profiles obtained from single evolution channels 
for comparison. 

Our results confirm the dependence of the pulsations on 
the shape of the C/O profile and indicate that binary evolu- 
tion may provide a better match to the pulsation frequencies 
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Fig. 4. The observed periods of GD 358 (dashed line), and the 
optimal model periods found by the genetic algorithm (solid 
line), using a profile for a remnant of 0.894 Mq from binary 
evolution with 2 CEs. 

of GD 358. Our best-fit model involves a previous evolutionary 
history with two CE episodes that has removed all of the hy- 
drogen from the surface of the star. This is an alternative way 
to explain the lack of hydrogen in DB white dwarfs. IT85 es- 
timated a formation rate of ~ lO'^yr'^ of non-DA degenerate 
dwarfs created from these binary scenarios, and concluded that 
they may produce 15%-20% of all non-DA white dwarfs. Our 
results suggests that GD 358 may be one of these. The best-fit 
model has a mass close to the mean mass for DB white dwarfs, 
a temperature consistent with recent calibration of UV spectra 
obtained with the lUE satellite, and a thin outer layer of helium. 

The value obtained for in our best-fit model (0.585 Mq) 
is lower than the IT85 mass value related to the profile that pro- 
duced this best-fit (0.894 Mq). In general, there are differences 
between these values for the different profiles tested. This may 
indicate that the physics of the evolutionary code and the white 
dwarf models are not completely consistent. However, in most 
of the cases, we have included profiles that use the same evo- 
lutionary code to be able to do qualitative comparisons, avoid- 
ing possible diff'erences in the fits due to the different physics 
applied to the models. Therefore, the method presented may 
also be used to investigate the coherence between evolution- 
ary codes and the adiabatic non-radial oscillation code. Further 
C/O profiles obtained with improved physics must be analyzed 
to investigate this relation. 

The results obtained for GD 358 indicate that this method 
may be a good option to investigate the evolution of DBV 
white dwarfs. This method is strongly related to observational 
properties of the stars, in particular to their pulsation periods. 
The main problem is that it can only be applied to few DBV 
white dwarfs: those that have been extensively observed in 
multisite campaigns in order to resolve their temporal spec- 
trum. Currently there are only two candidates suitable for this 
study: GD 358 and CBS 1 14 (Handler et al. 2002; Metcalfe et 
al. 2005). Further options for the future include extending this 
study using profiles coming from double evolutionary scenar- 
ios where the star experienced an epoch of mass loss from one 
CE episode during the AGB phase. 
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